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ABSTRACT 

We present diffraction-limited Ks band and L' adaptive optics images of the edge-on debris disk 
around the nearby F2 star HD 15115, obtained with a single 8.4 m primary mirror at the Large 
Binocular Telescope. At Ks band the disk is detected at signal-to-noise per resolution element (SNRE) 
~ 3-8 from ~ l-2'/5 (45-113 AU) on the western side, and from ~ 1.2-271 (63-90 AU) on the east. 
At V the disk is detected at SNRE ~ 2.5 from ~ 1-1745 (45-90 AU) on both sides, implying more 
symmetric disk structure at 3.8 /im. At both wavelengths the disk has a bow-like shape and is offset 
from the star to the north by a few AU. A surface brightness asymmetry exists between the two sides 
of the disk at Ks band, but not at L' . The surface brightness at Ks band declines inside 1" (~ 45 
AU), which may be indicative of a gap in the disk near 1". The Ks - L' disk color, after removal of 
the stellar color, is mostly grey for both sides of the disk. This suggests that scattered light is coming 
from large dust grains, with 3-10 /im-sized grains on the east side and 1-10 /mi dust grains on the 
west. This may suggest that the west side is composed of smaller dust grains than the east side, which 
would support the interpretation that the disk is being dynamically affected by interactions with the 
local interstellar medium. 

Subject headings: instrumentation: adaptive optics — techniques: high angular resolution — stars: 
individual (HD 15115) — circumstellar matter — planetary systems 
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1. INTRODUCTION 

Debris disks are thought to be signposts for planet 
formation since their reservoirs of dust grains must 
be frequently replenished by collisions between larger, 
planetesimal-sized bodies ()Wvattll2008D . The idea of de- 
bris disks as markers for planets has been supported 
by the recent direct detections of wide-orbiting, mas- 
sive exoplanets in systems that also harbor luminous de - 
bris disks fe.g..lMarois et all pOlOh : iKalas et"aTJ (|2008lh 
lLagrange et al.l (|2010)). In systems with spatially re- 
solved disks where no planets are detected, knowledge 
of dust grain sizes, compositions, density distributions, 
and disk morphologies can help us infer where any 
unseen planets mig ht reside (|Moro-Martfn fc Malhotral 
20051 [QuilleDl 12009) . 

To date the highest resolution images of debris disk 
scattered light have been obtained in the visible to near- 
infrared (NIR) with the Hubble Space Telescope (HST) 
and from the ground with adaptive optics (AO). Unfortu- 
nately, due to thermal emission from the Earth's atmo- 
sphere, obtaining ground-based detections of scattered 
light from extended debris disks at wavelengths longward 
of 2 fj,m is difficult. But the 2-5 [im wavelength region is a 
window into important debris disk properties. Scattered 
light at these wavelengths is sensitive to larger grains (> 
a few /um) than can be probed by visible to NIR scat- 
tered light. Imaging at these wavelengths can also con- 
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strain dust grain compositions, in particular the fraction 
of wa ter ice on the surface of dust grains ()Inoue et al.l 
2008). The 3-4 /xm wavelength region is also particu- 
larly sensitive to the therma l radia ti on of mass i ve exo - 
planets (e.g [Burrows et all (1200 3): Hi nz et all (|2010f ): 
IMarois et all (|2010h : iLagrange et al.1 ([20101)). Therefore 
imaging at these wavelengths simultaneously probes disk 
scattered light from large grains and thermal radiation 
from exoplanets. 

Because debris disks are optically thin, edge-on disks 
have a larger optical depth along the line of sight than 
face-on disks. Therefore, in the favorable case of an edge- 
on inclination, disk flux densities (assuming isotropic 
scattering) are maximized, facilitating their detection 
above the high sky background. In addition, when cre- 
ating a reference point spread function (PSF) from the 
intrinsic rotation of the sky relative to the telescope using 
the sc ience target (as in Liu (2004) and Fitz gerald et al.1 
(2007), for example), edge-on disks are more resistant to 
disk self-subtraction when differencing the reference PSF 
from the observed target. However being edge-on and 
brigh t does not gua rantee detections longward of 2 /xm, 
as iFitzgerald et al.l (2007) obtained confident detections 
of the AU Mic debris disk in the J, H, and Ks bands, but 
reached only a 3<r upper limit of ~ 12.6 mags/arcsecond 2 
at L' (in 12 minutes of integration) . To confidently detect 
scattered light at 3.8 /im from the ground, longer integra- 
tions (relative to integrations at shorter wavelengths) are 
required to mitigate the higher sky background. Addi- 
tionally, detections are facilitated if the disk consists of 
large, grey-scattering grains as opposed to small, blue- 
scattering grains. 

HP 15115 is a nearby (d = 45.2 ± 1.3 pc (jvan Leeuwenl 
120071 )) F2 star with an edge-on asymmetric debris disk 
previously spatial l y resolved i n the visible and NIR 
(jKalas et al.1 I2007t JDebes et aJ J 12008ft . The star is be- 
lieved to be young for several reasons: it has shared 
kinematics with t he 12 Myr old /3 Pic moving group 
(|Moor et al.l l2006h . it is believed to be on the zero-age 
main sequence (Eric Mamajek, private communication), 
and it has a high fractional luminosity circumstellar disk 
(fd = 4.9 x 1 Q~ ), which is m ore commonly seen for 
younger stars (M oor et "ail 120061) . However other indi- 
cators, such as Ca II H and K lines and X-ray emis- 
sion, may point to a much older age, p erhaps 100-500 
Mvr rtSilverstondlgOOCt iRhee et all 120071 ). Furthermore 
iDebes eta l (2008) refuted the evidence for the star be- 
ing a (3 Pic group member, based on backtracking the 
star's proper motion and radial velocity. Given the large 
uncertainty in this star's age, for the purposes of this pa- 
per we take the star's age to be conservatively between 
10 Myr and 1 Gyr. 

HD 15115's circumstellar disk is believed to 
be gas-poor, with < 0.28xl0~ 4 M e in CO gas 
(|Mo6r et al.ll2011aft. and has 0.047 M ffl in dust mas s 
(jZuckerman fc Song] 12001 I Williams fc Andrews! [2006h . 
Therefore the disk is considered to be predominantly 
"debris." At visible wavelengths, the disk is highly 
asymmetric, with the western lobe extending out to a 
stellocentric radius o f ~ 12" in a "need le" -like feature 
(jKalas et all 12007ft . iKalas et al.l (|2007l ) reported blue 
F606W - H colors, es pecially at la r ge sep arations (hence 
the "blue needle"), foebes et all (|2008| ) also saw blue 



scattering beyond 2", but reported red F110W - H colors 



at 1". The red color of the disk close to the star makes 
HD 15115 an attractive target for 2-4 /im imaging. 

Resolving the disk longward of 2 (im would allow us 
to constrain the population of large grains in the disk, 
especially closer to the star. This requires high-Strehl 
ratio, low thermal background observations. From the 
ground, we need a precise adaptive optics (AO) system 
and a minimal number of warm surfaces in the optical 
path. 

The Large Binocular Telescope (LBT) satisfies these 
requirements. Currently it has a single adaptive sec- 
ondary mirror with 675 actuators, capable of operating 
with up to 500 modes of correction on one of the two 
8.4 m primary mirrors. This allows for high-Strehl ratio 
(up to 70-80% at H band and up to 95 % at Ks band) 
low th ermal background observations (see Esposito et al. 
(2011) and references therein for a review of the LBT AO 
system) . 

We observed HD 15115 at Ks band and at L' with 
the DX (right) primary mirror and its adaptive sec- 
ondary mirror at the LBT. In Section 2 we describe 
the observati ons we carried out at Ks band with 
PISC ES dMcCarthvet alll200l and at V with LMIR- 
cam (Skrutsk ie et all [2.010) combined with the non- 
interferometric mode of the Large Bino cular Telescope 
Interferometer (LBTI. Irlinz et aLlpOO S)). as well as our 
data reduction methods. In Section 3 we present our re- 
sults on the disk structure and surface brightness (SB) 
profiles, and limits on planets in the system. In Section 4 
we discuss the implications of our results, including the 
disk color and grain sizes. In Section 5 we summarize 
the main results. 



2. OBSERVATIONS AND DATA REDUCTION 
2.1. LBT/PISCES Ks band 

We carried out our Ks band observations of HD 15115 
on UT November 9 2011 with the PISCES camera, a 
high-contrast 1-2.5 /an imager modified for use at the 
LBT. PISCES has a field of view (FOV) of - 19'.'5 on a 
side and a plate scale of 19.4 mas/pixel PISCES, which 
uses a pyramid wavefront sensor with natural starlight, 
was mounted at the right front bent Gregorian focus be- 
hind the First Light Adaptive Optics (FLAO) system 
(jEsposito et al.l l2010h . Skies were clear and the seeing 
was good (0.5-1") for most of the duration of the ob- 
servations. We used the single 8.4 m DX (right) mirror 
combined with its adaptive secondary mirror. The cam- 
era rotator was fixed at a fiducial position to allow f or 
angular differential imaging (API; IMarois et all J2006)), 
and no coronagraphs were used. We obtained 603 images 
with 4 s integration per image, resulting in a total inte- 
gration of 40.2 minutes and ~ 40° of parallactic angle 
rotation. For approximately the first half of the integra- 
tion, the star remained stationary on the detector and 
was not nodded or dithered. For the second half of the 
data, the star was moved to the opposite side of the de- 
tector, after approximately half an hour of down time 
during which the AO was not functioning. We also ob- 
tained several 0.8 s unsaturated images of HBC 388 in 
the narrowband BrG (Ao = 2.169 /im) filter, a few hours 
after HD 15115, as a photometric reference. 

The raw images were corrected for cross-talk and geo- 
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metrical distortiorQ . After these steps, all data reduc- 
tion was performed with custom Matlab scripts. Each 
half of the data was used to subtract the sky and detec- 
tor artifacts from the opposite half, and flat-field correc- 
tions were applied from sky images obtained earlier in 
the night. Each image was divided by its exposure time 
to obtain units of counts/s. All images were saturated at 
the cores (inside A/D) of the PSF, so the star positions 
were determined by calculating the center of light using 
an annulus from just outside the saturated region to 2". 
We performed similar center of light calculations on the 
unsaturated images of HBC 388 and found the difference 
between the true centroid and the calculated centroid to 
be ~ 1 pixel (0'.'194). The images were then registered 
to a fiducial position. For each image, we calculated the 
azimuthal median radial profile after masking out the 
spider arms and other detector artifacts, constructed an 
image built from the radial profile, and subtracted that 
from the original image. 

After median-combining the first half of the images, 
the second half of the images, and all images together, we 
determined that the second half of the data was of much 
poorer quality than the first half, perhaps due to the 
worse seeing during this part of the night. Consequently 
we only used the first half of the data in our final data 
reduction and analysis. We constructed the master PSF 
reference by median-combining the first 284 images, and 
scaling and subtracting this reference image from each 
of the first 284 images. The scale factor was determined 
by minimizing the flux after subtraction inside a 0.15- 
1" annulus (outside the saturated region). Each image 
was then flipped about the vertical axis (due to a mirror 
flip inside PISCES) and rotated by its parallactic angle 
plus an additional rotational offset to obtain North-up, 
East-left. The offset was determined to be 80.88° (the 
position angle of the detector) plus 3.9° (see Section 3.2 
for a discussion of how the 3.9° offset was determined). 
The total parallactic angle rotation for this set of images 
was approximately 20°. 

We median-combined the PSF-subtracted, rotated im- 
ages, masking out the residual streaks left by the spider 
arms and other detector artifacts. The final image re- 
vealed HD 15115's edge-on asymmetric disk structure at 
the expected posit i on an gle (PA) of ~ 279° (Kal as et al.l 
[20071 iDebes et al.l l200g )). To obtain the highest SNR 
image, we reduced the first 284 (better quality) i mages 
using a cons ervative LO CI (Lafrcnierc et al. 2007) algo- 
rithm , as in IBuenzli et al.l (|2010f ) and iThalmann et al.l 
(|2011l) . Specifically, to preserve disk flux we required 
a field rotation of at least 3 times the full-width half- 
max (FWHM) at Ks band (= 3 X 0'.'0525) between im- 
ages. We also set the optimization section (Na) to 2000 
to ensure less self-subtraction. Fig. [Ta| shows our final 
PISCES/As band image in mJy/arcsecond 2 . 

Because we could not obtain any unsaturated images 
of HD 15115 (the shortest integration time possible with 
PISCES is 0.8 s, which saturates the star), the flux cal- 
ibration for obtaining units of mJy/arcsecond 2 was cal- 
culated as follows: we scaled the ratio of the PSF halo 
of HBC 388 to its total flux and compared that to the 
halo of the median-combined PSF reference image of HD 
15115 (so that the disk flux contribution was washed 

1 see http://aries.as.arizona.edu/ and links therein 



out). The ratio was determined for different halo sizes, 
ranging from 0'.'3 to 1", l'/25, and 1'/ 5— all yielding the 
same ratio to within less than 1%. We adopted this 
method, as opposed to scaling by the relative bandwidths 
in the Ks and BrG filters, because the latter method 
may be less accurate due to changing photometric con- 
ditions. 

Fig. Ilblshows the map of the signal-to-noise per reso- 
lution element (SNRE) of the final Ks band image. The 
SNRE was determined by masking out the disk in the 
final image, smoothing the image by a Gaussian kernel 
with FWHM = 0'.'0525 (= A/D), calculating the stan- 
dard deviation of the smoothed image in concentric an- 
nuli around the star, then dividing this noise image into 
the final smoothed Ks band image. From the SNRE 
map, the disk is detected on the eastern side at SNRE 
~ 3 between 1.2-2'.'l; and on the western side the disk is 
detected at SNRE ~ 3-8 between l-2'/5. 

We refer the reader to the Appendix, wherein we de- 
scribe how we inserted an artificial model disk into the 
raw data and re-reduced the data using the LOCI al- 
gorithm, to understand LOCFs effects on disk position 
angle, FWHM, and surface brightness as a function of 
distance from the star. The algorithm's effects on these 
values were measured and accounted for in the disk anal- 
ysis found later in this paper. 

2.2. LBTI/LMIRcam V 

We carried out our V observations of HD 15115 on 
UT N ovember 16 2011 with LMIRcam (jSkrutskie et al.l 
120101) . a hig h-contrast 3-5 jiy a. imager designed for use 
with LBTI (jHinz et al.ll2008[ ). These observations were 
made in direct imaging mode, without any interferomet- 
ric combination. Like the FLAO system, LBTI also uses 
a pyramid wavefront sensor with natural starlight for 
the adaptive reference source. We used the 8.4 m DX 
mirror combined with its adaptive secondary. LMIR- 
cam has a FOV of ~ 11" on a side and a plate scale of 
10.7 mas/pixel, determined from observations of the bi- 
nary star HD 37013 the previous night. LBTI is fixed 
to the telescope and thus cannot rotate, facilitating the 
use of ADI throughout the observations. We obtained 
2319 images, the first 5 with integrations of ~ 10 sec- 
onds, and the remainder with integrations of ~ 2.5 sec- 
onds. The total integration obtained on the source was 
therefore ~ 97 minutes. Observing conditions were good, 
with clear skies and ~ 1" seeing. For much of the 
second half of the integration, the wind picked up and 
caused some vibration-induced blurring of the saved im- 
ages. At the end of the acquisition, we obtained several 
short-exposure, unsaturated images of H D 1511 5 and HD 
40335, an L' standard star (Lcggctt et al. 2003), for pho- 
tometric references. Throughout the long exposure im- 
ages, each image was saturated inside A/D. Due to a 
poorly-mounted dichroic, astigmatism in the PSF was 
visible inside the central ~ 0'.'2 of each image, creating a 
"cross" pattern on the PSF. This did not ultimately have 
any measurable effect on the quality of the data taken, 
and the instrumental astigmatism was corrected after the 
observing run. Throughout all observations, we nodded 
the star diagonally by several arcseconds, and dithered 
the star by 0'.'25 in a box around the center of each nod 
position. 

All data reduction was performed in custom Matlab 
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Fig. 1.— Top: Final PISCES #"s band image of the HD 15115 
debris disk, in units of mjy/arcsecond 2 , with North-up, East-left. 
The white dot marks the location of the star, and represents the 
size of a resolution element at Ks band. For display purposes a 0'.'8 
radius mask has been added in post-processing. The western side 
SB is ~ a magnitude/arcsecond 2 brig hter than the eastern side SB, 
as is seen at shorter wavelengths (Kalas ct al. 2007; Dcbcs et al. 
120081) . Bottom: SNRE map of the image. The eastern side is 
detected at SNRE ~ 3 out to ~ 2'.' 1 . The western side of the disk 
is detected at SNRE ~ 3-8 from ~ l-2'.'5. 

scripts. Images were divided by their individual integra- 
tion times to obtain units of counts/s, nod-subtracted, 
and corrected for bad pixels. Due to the large quantity 
of data, images were averaged in sets of 4. We registered 
the images to a fiducial pixel by calculating the center 
of light around the star, excluding the saturated pixels 
within ~ A/D, in the same manner as for the Ks band 
images. We performed similar center of light calcula- 
tions on the unsaturated images of HD 15115 and found 
the difference between the measured and true centroid to 
be ~ 0.2 pixels (2 mas). We then calculated the radial 
profile of each image and subtracted it, also in the same 
manner as for the Ks band images. To test the quality of 
the data, we median-combined the first half of the data, 
the second half, and the entire set. We saw that though 
the quality of the PSFs was good, the second half of the 
data suffered from bright extended streaks appearing in 
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Fig. 2. — Top: Final V image of the disk obtained with LMIR- 
cam, in units of m Jy/arcsecond 2 , with North-up, East-left. The 
white dot marks the location of the star, and represents the size of 
a resolution element at L' . For display purposes a 0'.'9 radius mask 
has been added in post-processing. The image has been smoothed 
by a Gaussian kernel with FWHM = A/D and binned by a factor 
of 4 to bring out the disk. At L' the disk is mostly symmetric and 
is equally bright on both sides. These features are nearly opposite 
to what is seen at Ks band and at shorter wavelengths. Bottom: 
SNRE map of the final L' image shown above. The disk is detected 
at SNRE ~ 2.5 between ~ l-T'5 on both sides. 

the raw data, out to several arcseconds. These streaks 
were not caused by the spider arms, and were not sym- 
metric, so they could not be easily masked out and would 
overwhelm the brightness of any disk structure seen in 
the final image. Therefore we used only the first half 
of the data (~ 40 minutes of integration, with ~ 50° of 
parallactic angle rotation). 

After scaling and subtracting (in the same manner as 
for the Ks band images) the master PSF from each im- 
age in the first half of the data (240 averaged images; 
960 total), we rotated each image by its parallactic an- 
gle plus an offset to obtain North-up, East-left. The 
offset was determined to be 1.81° ± 0.0685° by calculat- 
ing the PA of the binary HD 37013, which was observed 
with LMIRcam the night before. The error was deter- 
mined by independently calculating the offset with two 
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different software routines and measuring the difference. 
The PSF-subtracted images were then median-combined. 
The final image, after smoothing by a Gaussian kernel 
with FWHM = A/D and binning by a factor of 4, re- 
vealed edge-on disk structure at a similar PA to the disk 
in the PISCES Ks band image. 

To obtain the highest-SNR image, we reduced the data 
using a conservative LOCI algorithm. Here we required 
3 FWHMs (= 3 x O'.'l) of rotation between images, and 
the optimization section Na was 300 (~ l'/5). Fig. [2a] 
shows this final image in mJy/arcsecond 2 , where the con- 
version to mJy/arcsecond 2 was dete rmin ed using the un- 
saturated images of HD 15115. Fig. l2bl shows the SNRE 
map of the final L 1 image, where the SNRE map was cal- 
culated in the same manner as for the Ks band image. 
The disk at V is detected at SNRE ~ 2.5 from l-T/5 on 
both sides of the disk. 

We refer the reader to the Appendix, wherein we de- 
scribe how we inserted an artificial model disk into the 
raw data and re-reduced the data using the LOCI al- 
gorithm, to understand LOCFs effects on disk position 
angle, FWHM, and surface brightness as a function of 
distance from the star. The algorithm's effects on these 
values were measured and accounted for in the disk anal- 
ysis found later in this paper. 

2.3. Ancillary Keck/NIRC2 Ks band data 

During the LBT observations, the PISCES camera 
lacked an accurate astrometric solution to determine the 
PA of true North on the detector. This is largely due to 
faulty fits header reporting of the detector PA, as well as 
the lack of observations of a standard binary for astro- 
metric calibration. Therefore we compared our PISCES 
final image to Ks band data of HD 15115 obtained on the 
night of UT August 12 2011 with Keck/NIRC2, which 
has a very precise astrometric calibration and distortion 
solution determine d from repeat ed observations of the 
galactic center (fYelda et a.l.ll2010h . 

Briefly, the data were taken with the narrow camera 
with the 0'.'6 diameter coronagraphic mask and consisted 
of coadded 10 s exposures, with a cumulative integra- 
tion time of ~ 1280 seconds and total field rotation of 
~ 39°. The data were processed using a conservative 
LOCI algorithm. To provide a detection of the HD 15115 
disk minimally-biased by LOCI processing, we adopted 
relaxed LOCI paramet ers comparable to those used in 
iThalmann et al.l (|2011[ ). Specifically, we set the LOCI 
parameter for the required rotation gap, 5, to be 3 times 
the FWHM (3 x 070 52) and set the opt imization area, 
N A , to be 3000 (see iCurrie et al.1 (J2011D for additional 
details on the algorithm used). The disk was detected 
between ~ 1-2" at SNRE ~ 3-5. The final image was ro- 
tated clockwise by an additional 0.252° to obtain North- 
up, East-left. We show only the PISCES Ks band image 
in this paper because the disk in the PISCES image was 
detected at higher SNRE. The Keck/NIRC2 data will be 
discussed and analyzed in Currie et al. (2012, in prep.). 

3. RESULTS 

3.1. Disk FWHM 

To better define the Ks band morphology of the disk, 
we measured the disk FWHM perpendicular to the major 
axis of the disk as a function of stellocentric distance. 



We did this as follows: first we rotated the final image 
by 7° so that the disk midplane was horizontal; next at 
each discrete horizontal (pixel) distance from the star 
we located the brightest pixel in the column of pixels 
perpendicular to the disk; we then placed a 7 pixel by 
25 pixel (0'.'13 by 0748) box centered on the brightest 
pixel, and summed up the counts/pixel along each row 
of the box, producing a ID array of 25 values. Next 
we fit a Gaussian to this array. The Gaussian fit gives 
the disk midplane pixel location along with the FWHM 
at that location. We also measured the disk FWHM 
as the number of pixels along each row with SNRE > 
3. We took the final disk FWHM value as the average 
of the Gaussian-fitted value and the SNRE- width value, 
with the error being half the difference between the two. 
We also added in the constant FWHM offset (~ 0703) 
determined from insertion and recovery of an artificial 
disk in the raw da ta (s ee Appendix for details). 

Fig. I3"al and Fig. I3bl show the east and west FWHM of 
the disk as a function of distance from the star, respec- 
tively. The western disk FWHM increases gradually with 
increasing distance from the star, though this may be an 
effect from the LOCI reduction (see Appendix), while no 
strong correlation is evident for the eastern FWHM. The 
median FWHM across the entire disk was found to be 
0721 ± 0703. The intrinsic FWHM of the disk is com- 
puted by subtracting off in quadrature the PSF broaden- 
ing at Ks band (= 070525), but this has a negligible ef- 
fect on the disk FWHM, yielding the same value. Within 
the uncertainty, this value is consistent with both the in- 
trinsic disk FWHM at 1.1 /im, 0726 ± 0'.'07, where we 
computed this value by averaging th e reported intrinsi c 
disk FWHM of each side of the disk (|Debes et al.ll2008j) . 
as well as the FWH M of the disk at 0.6 /im, reported by 
iKalas et al.l (|2007l ) as 0719 ± O'.'l. The measured disk 
FWHM at Ks band being > 3 times the size of the PSF 
FWHM (070525) suggests that we have spatially resolved 
the disk. 

We also measured the FWHM of the disk at L'. We 
did this in the same manner as for the Ks band image, 
except the L' image was binned by a factor of 4 to in- 
crease the disk SNR per pixel, and we used 3 pixel by 
11 pixel (~ 0713 by 0747) boxes. We also added in the 
constant FWHM offset (~ O'.'ll) determined from inser- 
tion and recovery of an artificial disk in the raw data 
(see Appendix for details). Given the low SNRE (~ 2-3) 
of the L' disk and the lack of spatial coverage (detected 
between only 1-175), we do not plot the FWHM as a 
function of distance from the star. Instead we computed 
the global FWHM value as the median of the measured 
FWHM disk values, which was found to be 0726 ± 0703. 
The intrinsic disk FWHM, calculated by subtracting off 
in quadrature the PSF broadening at L' (= 070940), is 
0724 ± 0703. Within the uncertainty this value is consis- 
tent with the both the intrinsic disk F WHM at Ks band , 
as well as the di sk FWHM at 1.1 /xm (]Debes et al.ll2008h 
and at 0.6 /im (|Kalas et al. 2007]). Given that the mea- 
sured FWHM is ~ 2 times the size of the PSF FWHM 
at L', the disk is marginally spatially resolved. 

3.2. Disk Position Angle 

To determine the rotational offset needed to obtain 
North-up, East-left in our final PISCES Ks band im- 
age, we compared our Ks band disk image to the 
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Fig. 3. — Top: Ks band eastern disk FWHM as a function of 
distance from the star. See text for methodology. There is no 
obvious trend in FWHM vs. stellocentric distance. Bottom: Same 
as above, except for the western side of the disk. The FWHM 
increases with increasing distance from the star, though this may 
be an effect from the data reduction (see Appendix). 

Keck/NIRC2 Ks band image, since the rotational off- 
set has been well-calibrated for Keck/NIRC2. We do not 
compare our PISCES image to published HST images 
to avoid any potential wavelength-dependent changes in 
PA. Comparing our final PISCES Ks band image with 
the Keck image between 1-2" revealed a difference in 
overall PA of 3.9° , so we rotated our PISCES image by an 
additional 3.9° clockwise. Fig. [la] shows this correctly- 
rotated image. 

To test for any warps in the disk, we calculated the PA 
of the disk at Ks band on both sides of the disk as a func- 
tion of distance from the star, using the disk midplane 
pixel locations from the Gaussian fits. The errors were 
measured as the difference between these values and the 
center of light along each ID row at each discrete pixel 
distance from the star. The errors also included the 1 
pixel = 0'.'194 astrometric centroiding uncertainty. Fig. 
[4al and Fig. [4b] show the Ks band PA of the eastern and 
western sides of the disk as a function of distance from 
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Fig. 4. — Top: PA of the eastern side of the disk at Ks band 
as a function of distance from the star. See text for methodology. 
The PA increases with increasing distance from the star; over the 
same stellocentric distances a similar trend was observed at 1.1 
/mi (Dcbcs ct al. 2008). Bottom: Same as above, except for the 
western side . The PA increases closer to the star, as is seen in the 
IDebes et all J2008I) 1.1 /im band data. 

the star, respectively. The eastern PA increases with in- 
creasing distance from the star, similar to what is seen 
over the same spatial locations at 1.1 ^m ([Debes et al.l 
2008). The western PA increases closer to the star. This 
also agrees with the western PA vs. stellocentric distance 
seen at 1.1 ^m (|Debes et al.l 120081 ) and with the overall 
PA (278.5°) measured fo r the western side of the disk at 
0.6 /im between 1.5-12" (jKalas et alJl2007| ). 

We calculated the PA of the disk at V in the same 
manner as for the Ks band image, except we used 3 
pixel by 11 pixel (~ 0"13 by 0"47) boxes. The V image 
was again binned by a factor of 4. We do not plot PA vs. 
distance from the star due to the poor spatial coverage 
of the disk and, even after binning, the low SNRE (~ 2- 
3); however the measured values for both sides of the 
disk were globally consistent with the values measured 
at Ks band (this is also evident by examination of the 
Ks band and L' disk images in Fig. [la] and Fig. [2aJ 
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respectively) . 

3.3. Surface Brightness Profiles 

For both the Ks band and the L' disk images, we ob- 
tained the SB profiles by calculating the median disk flux 
in circular apertures with a radius of 0'.'15 centered on 
the Gaussian-fitted disk midplane location at each dis- 
crete pixel distance from the star. We chose this aperture 
size because it maximized the SNR of the disk flux at L'. 
This was necessary because the L' disk is detected at 
low SNRE (~ 2.5), so using larger apertures incorpo- 
rated more background noise. We measured the errors 
as the standard deviation of the median SB values in the 
same apertures placed in a circular annulus at each pixel 
distance from the star. We computed the aperture cor- 
rection at Ks band to be 0.64 by convolving the image of 
the photometric reference HBC 388 with a bar of unity 
counts/pixel and ~ equal width to the observed disk at 
Ks band. The SB as a function of distance from the star 
is shown in Fig. [5a] "Magnitudes" in this paper refers 
to Vega magnitudes. 

Between ~ 1.2-2'.'l (where both sides of the disk 
are detected), the typical western disk SB is ~ a 
magnitude/arcsecond 2 brighter than the eastern side. 
A s imilar bri ghtness asymmetry was seen at 1.1 
jum (|Debes et a l. 2008). It is also interesting that the 
western SB appears to drop off near 1", by ~ half a 
magnitude/arcsecond 2 (although the noise is high inte- 
rior to this distance). We saw this SB reduction in our 
final ADI image as well as in our final LOCI image (Fig. 
Ilap . even after correction with insertion and recovery of 
an artificial disk. The d rop in SB near 1" was not seen 
at sh orter wavelengths (jKalas e t al. 2007; Dcbe s et al.l 
2008), which raises suspicion that the feature may not 
be real. However the shorter wavelength results appear 
to be limited by PSF residuals interior to 1", so it may 
not be appropriate to compare to these data. 



TABLE 1 

Ks BAND SB POWER-LAW INDICES 



Side 



index 



West (l-l'/2) 
West (1'.'2 < r < l'/8) 

West (r > l'/8) 
East (1'.'2 < r < 1'.'4) 

East (r > l'/4) 



3.48 ±0.5 
-1.75 ±0.05 
-4.40 ±0.24 
-5.27 ±0.41 
-1.1 ±0.07 



Table Q] shows the power- law indices measured for the 
eastern and western sides of the disk at Ks band. The 
power-law indices do not agre e very well with the in- 
dices at 1.1 /im reported by iDebes et al.1 ((2008). The 
discrepancies could be explained by the differing spatial 
coverage of the disk at the two wavelengths. On the 
other hand, it is also possible that the disk SB falls dif- 
ferently at Ks band than at 1.1 /im. Additional high- 
contrast imaging data at these wavelengths would help 
clarify this matter. We do not calculate the power-law 
indices at V due to the lower SNRE and poor spatial 
coverage of the disk. 

Using the Ks band power-laws, we can quantify the 
significance of the SB reduction interior to l'.'l. From 
the power-law index for l'/2 < r < l'/8, the expected 
disk SB should be 13.48 mags/arcsecond 2 at ~ 1". The 
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Fig. 5. — Top: SB profile of the disk at Ks band. The western 
side of the disk is ~ a magnitude/arcsecond 2 brighter than the 
eastern side at Ks band over the same stclloccntric distances. The 
western SB profile shows a ~ 2<r drop off in SB near 1-1" 2. This 
cannot simply be explained by flux loss due to disk self-subtraction, 
since this has been corrected for. Bottom: SB profile of the disk 
at L' . Both sides of the disk are ~ equally bright beyond 1'.'2, 
and just as in the Ks band image there is low SNR evidence for a 
decrease/flattening in SB near l-l'/2. 

actual value is ~ 13.82, with a la upper limit of 13.64. 
This suggests that we are seeing a ~ 2a reduction in SB 
interior to l'.'l. 

We calculated the SB vs. distance from the star for the 
L' image in the same manner as for the Ks band image, 
using the same size aperture as we did for the Ks band 
image. Additionally the image was binned by a factor 4 
to increase the signal-to-noise per pixel. We computed 
the aperture correction at L' to be 0.63 by convolving 
the photometric image of HD 15115 with a bar of unity 
counts/pixel and ~ equal width to the measured disk at 
V . We computed the errors in the same manner as for 
the Ks band SB profile. Fig. I5bl shows the SB profile of 
the disk at L' . Both sides of the disk are ~ equal in SB 
beyond 1'.'2, which is not seen at Ks band. Interestingly, 
within the uncertainties the data are consistent with a 
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flattening in SB interior to ~ l'.'l, similar to the SB drop 
off interior to l'.'l seen in the Ks band image. 

3.4. Limits on Planets 

We do not detect any high SNR point-sources in any 
of our data reductions that would point to a possible 
planet. To set firm constraints on what we could have 
detected, we reduced the same data that was used to de- 
tect the disk at L' with an aggressive LOCI algorithm. 
We examined only the V data, as opposed to includ- 
ing the Ks band data as well, because we had nearly 
double the integration time at V (~ 40 minutes vs. 20 
minutes at Ks band) and double the parallactic angle 
rotation (~ 40° vs. ~ 20° at Ks band). For the ag- 
gressive LOCI reduction we required only 0.75 FWHM 
of parallactic angle rotation between images, and our op- 
timization section size (-/Va) was 100. Our final image (in 
units of <7, where a was calculated in concentric annuli 
the size of the PSF FWHM (~ 0."1)), shown in Fig. El 
reveals no 5a point-source detections. 

To ascertain the minimum flux density point-source ob- 
ject we could have detected, we inserted artificial planets 
into the raw data and re-ran our aggressive LOCI reduc- 
tion. The artificial planets were bright PSFs obtained 
from the unsaturated images of HD 15115. We did not 
scale down the PSFs in brightness. Instead we simply 
calculated the SNR for each artificial planet's location 
(inside a radius = FWHM aperture) and determined the 
flux that would ha ve r esulted in a 5a detection. These 
are plotted in Fig. I6bl We reach a 5a background limit 
of ~ 17.8 mags (in 40 mins of integration). At 0'.'6, the 
contrast (computed by subtracting the star's L' magni- 
tude of 5.763) is ~ 9.5 mags. These sensitivity and con- 
trast limits show the incredible potential of direct imag- 
ing with LBTI/LMIRcam. For comparison, with Clio-2 
at the 6.5 m MMT, a 5a background limit of 17.8 mags at 
L' wa s achieved in 2.5 hours of integration (|Rodigas et al.l 

unni). 

We also overplot the magnitudes of 1-3 Mj plan- 
ets (10 Myr old) and, conservatively, the magni- 
tud e of a 30 Mr br own dwarf (1 Gyr old) from 
the lBaraffe et a l. (2003) COND mass- luminosity models 
(horizontal dashed lines). Assuming a young stellar age, 
at 5a confidence we rule out planets more massive than 
~ 1-2 Mj outside of 1", and outside of 0'.'6 we can rule 
out planets more massive than ~ 3 Mj. If HD 15115's 
age is much older, then we can only rule out brown dwarfs 
more massive than ~ 30 Mjbeyond 0'.'7. 

4. INTERPRETATIONS 

4.1. Disk Structure 

It has been suggested that HD 15115 is interacting 
with the local interstellar medium (ISM) , given its space 
motion to the s outh-east nearly p arallel with its disk 
major axis PA (jDebes et al.l I2009T ) a nd the disk asym- 
metries seen at sh orter wavelengths (|Kalas et al.l 120071 : 
iDebes et al.l [2008). The ~ 1 magnitude/arcsecond 2 
brightness asymmetry (between 1.2-2'.'l) and the east- 
west morphological asymmetry seen in our Ks band data 
both support this proposition. These effects could be 
caused by the eastern side of the disk plowing head-first 
into the ISM. In this case the eastern side would be much 
more affected than the western side, and small grains 
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Fig. 6. — Top: Final image at L' , in units of cr, computed with an 
aggressive LOCI reduction. See text for methodology. The white 
dot marks the location of the star and the size of a resolution ele- 
ment at L' . No high SNR point-sources are detected. Bottom: 5a 
sensitivity curve for our L' data, computed by inserting artificial 
planets into the raw data and running our aggressive LOCI reduc- 
tion. We reach a background limit of ~ 17.8 mags, and a contrast 
of ~ 9.5 mags at ~ 0"6. The horizontal dashed lines represent 
the magnitudes of 1 Myr old and 10 Gyr old planets, from the 
Baraffc et al. (2003) COND mass-luminosity models. We rule out 
planets more massive than ~ 3 Mj if the star is young and 30 
M j if the star is old. 

could be blown out to the west. This could cause the 
eastern side to be more truncated and fainter relative to 
the western side, both of which we observe. Any large (> 
a few /im) grains in bound orbits should remain mostly 
unaffected. Since the disk is symmetric at V , which is 
sensitive to large, grey-scattering grains, our data sup- 
port the ISM-interaction interpretation. 

It is also certainly possible that other dynamical ef- 
fects (e.g., pl anets or a close pa ssage of a nearby star as 
suggested bv iKalas et al.l (|2007l )) are responsible for the 
observed morphological and SB asymmetries. However 
there is currently no supporting evidence for a stellar 
flvbv (IKalas et al.1 120071 [Debes et all 12051 12001 . Dy- 
namical modeling of the disk with embedded planets is 
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Fig. 7. — Top: Rotated model and real Ks band disk images. 
The model image has an inclination of 87°, scattering parameter 
g = 0.5, 3 fiva grains, a SB power-law index of -3/2, and a gap from 
0-1'.' 1. The model image has been convolved with the PSF from 
the photometric image of HBC 388, and Gaussian noise has been 
added. Units are in mjy/arcsecond for both images. The white 
dot marks the location of the star in both images. The dashed lines 
are meant to guide the eye to the apparent northern offset of the 
disk relative to the star. Bottom: PA vs. stelloccntric distance for 
the western side of the disk in the model image. The PA increases 
closer to the star due to the apparent bow -shape, similar to what 
is seen for the real data shown in Fig. I4bl 

beyond the scope of this paper, but such modeling would 
help clarify what exactly is causing the observed asym- 
metries in the disk. 

Inspection of our Ks band and L' disk images (Fig. Hal 
and l2"a)) . as well as the PA vs. stellocentric distance plots 
for both sides of the disk (Fig. Hal and lib! suggests that 
the debris disk has a bow-like shape and is offset by a few 
AU to the north from the star. Specifically the offsets 
are ~ 0"012 (0.5 AU) and 0'.'04 (2 AU) at L', for the east- 
ern and western sides, respectively; and ~ O'.'ll (5 AU) 
and 0'.'17 (8 AU) at Ks band for the eastern and western 
sides, respectively. Since the astrometric uncertainty in 
centroiding is ~ 0'.'002 and 0'.'0194 at L'and Ks band, re- 
spectively, it is unlikely that the northern offsets could be 



due to centroid error. Bow-like disk morphologies have 
been seen in systems moving in near-perpendicular direc- 
tions to their disk major axis positi ons (e.g., HD 61005 
(jHines et all 12007ft and HD 32297 (|Debeset all [2009)1 
Because HD 15115 is moving to the south-east, nearly 
parallel with its disk major axis, one would not expect 
the disk to have a bow-shape, especially with an offset 
from the star to t he north. How e ver, in the HD 61005 de- 
bris disk system, iManess et all (|2009|) showed that dust 
grain interactions with the ISM in directions along the 
disk midplane could actually create bow-shapes perpen- 
dicular to the disk major axis, roughly reproducing that 
disk's observed morphology. Similar dynamical interac- 
tions in the HD 15115 system could be creating the ob- 
served bow-shape, though no explicit model has tested 
this hypothesis for this system. 

Another perhaps simpler explanation for the bow- 
shape and the apparent offset is that these are caused by 
geometrical effects. We can reproduce the observed bow- 
shape and vertical offset with a simple inclined, ringed 
disk model. We used the radiative transfer equations 
describing the intensity of light scattering off dust par- 
ticles and made three assumptions about the disk: it is 
not exactly edge-on, but rather inclined to 87°; it has 
a gap from 0-1'.' 1 (see Section 4.3 for additional discus- 
sion of the possible disk gap); and we set the Henyey- 
Greenstein scattering parameter g = 0.5 for predomi- 
nantly forward-scattering grains. This is a reasonable as- 
sumption for large (> 1 /xm) grains scattering NIR light 
(jDuchene et all 120041 ) . We set the dust grain size uni- 
formly to 3 (im and set the disk SB power-law index to 
-3/2. Fig. ITal shows the rotated model image, along with 
the rotated final Ks band image. Dashed lines have been 
inserted to guide the eye to the app are nt bow-shape and 
northern offset from the star. Fig. [7b] shows the PA vs. 
stellocentric distance for the western side of the disk in 
the model image, computed in the same manner as for 
the real data. The PA increases closer to the star, as is 
observed in the real data (Fig. I4b|) . 

While the parameters used in our model are not a 
unique explanation for the observed features of the disk, 
the proposition that the effects can be explained by the 
disk's geometric orientation in space is attractive because 
it explains the observations without contradicting the 
evidence supportin g the ISM interaction interpretation 
(|Debeset al.ll2009T ). 

4.2. Disk Color and Grain Sizes 

To determine the disk color as a function of distance 
from the star, we calculated A(Ks - L 1 ) — (Ks - L')nisk 
- (Ks - L')star m tie regions of spatial overlap between 
the Ks band and V images (1.1-1'.'45). The errors were 
calculated by summing the individual Ks and L' errors 
in quadrature. The data suggest that the eastern side 
of the disk becomes redder than the western side with 
increasing distance from the star. Interior to 1'.'3 both 
sides of the disk are grey. 

To constrain the characteristic dust grain sizes, we 
compared our A(K s - L 1 ) disk colors to the models of 
llnoue et"a l (2008). These models calculate disk colors 
in the NIR and assume silicate dust grain composition. 
The modeled grain sizes range from 0.1 //m to 10 /im and 
assume a single grain size for the dust population. We 
plot the model disk colors for 1, 3, and 10 fim grains. 
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Fig. 8. — Disk color vs. distance from the star, expressed as 
A(Ks - L') = (Ks - L')rj)i s fc - (Ks - L')star- This is calculated only 
where the disk is detected at both Ks band and at L'. To constrain 
dust grain size, we also plot model colors fromjlnouc ct al] 1120081 ) 
(colored horizontal lines). The data suggest that the eastern side 
of the disk becomes redder than the western side with increasing 
distance from the star. 1-10 fim grains are the best fit to the 
western side, while 3-10 /an grains are the best match to eastern 
side. This may suggest that the west side of the disk is composed 
of smaller grains than the east side. 

These are shown as the colored lines in Fig. [51 We also 
modeled 0.1 /im and 0.3 /im grains, but their colors were 
too blue to be supported by the data. 

From Fig. we see that 1-10 /mi-sized grains are the 
best match to the western side of the disk, while 3-10 
/itm are the best match to the eastern side. This sug- 
gests that the observed scattered light from the disk is 
tracing the large parent body dust grains in the disk. 
The data also may suggest that the west side of the disk 
is composed of smaller grains than the east side. The 
blowo ut grain size for this system is expected to be ~ 1- 
3 fim (Ha hal2010l ). assuming a stellar mass of 1.3 M©, 
a luminosity of 3.3 L©, and a grain density in the range 
1-2.5 g/cm 3 . Grains smaller than the blowout size would 
be blown out radially; on the eastern side, these grains 
would hit the ISM and be blown back to the west, re- 
sulting in the western side being dominated by smaller 
blue-scattering grains. Our observational constraints on 
the dust grain sizes offer some support for these predic- 
tions. 

4.3. Does the disk have a gap? 

The SB profiles at Ks band and at L' drop off or flatten 
out (at low SNR) near 1". For an edge-on disk with no 
gap, the SB should continue to increase closer to the 
star. Because we do not see this in our data, this may be 
an indication that the disk has a gap interior to ~ l'.'l. 
This is con s istent w ith prediction of a gap near 1" by 
IMoor et all (|2011bl ). using spectral energy distribution 
(SED) analysis of the system. Though degenerate with 
temperature and dust grain size, their best-fit model of 
HD 15115's infrared SED yields a two-component disk, 
with the inner at 4 AU and the outer at 42 AU (both ±2 
AU). Given a distance to the star of 45.2 pc, the outer 
edge of the gap would be at 0'.'93 ± 0'.'04. This agrees 
well with the observed drop offs in SB near 1". 

We can independently constrain the parent body dust 
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Fig. 9. — Equilibrium disk temperature vs. distance from the 
star, for several different silicate grain sizes. See text for methodol- 
ogy. The horizontal dash ed lines r epresent temperature constraints 
of 179 K and 57 K from IMoor et al.l l |2011bl ). Given that we have 
observational evidence for the predicted gap at ~ 1", we can inde- 
pendently constrain the dust grain sizes to be ~ 3 fim for astro- 
nomical silicates. 

grain size using estimates of the disk's equilibrium tem- 
perature as a function of distance from the star, for differ- 
ent grain sizes (Fig. [9]). The equilibrium dust tempera- 
ture was computed by balancing the absorption and the 
emission energy of a particle wit h adopted dust prop - 
erties for astronomical silicates (jLaor fc Draind 11993). 
Under an optically thin condition, the heating is solely 
from the central star, representative as the best-fit Ku- 
rucz model with a stellar temperature of 7000 K and a 
stellar luminosity of 3.3 L©. 

The horizontal dashed lines represent th e predicted 
disk tem peratures of 179 K and 57 K from IMoor et al.l 
(20 lib]). Assuming we have observationally detected the 
inner edge of the outer disk component (and hence the 
gap) at ~ l'.'l, the parent body dust grain sizes are con- 
strained to be ~ 3 ^tm. This is consistent with our dust 
grain size estimates from the disk colors. 

4.4. Limits on a planet inside the gap 

We can also calculate an independent estimate of the 
mass of a planet creating the gap in the disk using the 
equation describing the relationship between the width 
of the chaotic zone around an assumed coplanar, low- 
eccentricity pla net and its mass and semimajor axis 



eccentric ity pla ne 
(|Malhotralll998h : 



Aa w 1.4 a p (M p /M*f'\ 



(1) 



where Aa is the width of the chaotic zone, a p and M p are 
the semimajor axis and mass of the planet, respectively, 
and M* is the mass of the star. We make the assumption 
that, from our SB profiles, the outer edge of the gap, 
r gap , is at l'.'l (50 AU). Since Aa = r gap — a p , we can 
substitute this into Eq. [TJand solve for the planet mass 
as a function of semimajor axis (Fig. ITU]) . 

The planet's mass approaches zero the closer it orbits 
to the disk edge and approaches the brown dwarf mass 
regime closer to the star. Assuming the gap was created 
by an object whose semimajor axis = its projected sep- 
aration at the epoch of our observations and that the 
system is young, we can say that the object must be < a 
few Mj because we could have detected > 3 Mj planets 
(dashed line in Fig. Ilpp . The planet's semimajor axis 
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ffm (Debes et al. 2008) and at 0.6 /im (|Kalas et al.1 




BiiirT i 



30 M (1 Gyr) observational constraint 




semimajor axis (AU) 

Fig. 10. — Masses of a possible object orbiting inside a disk gap 
with its outer edge at l'.'l (50 AU), as a function of semimajor axis, 
computed using Eq. [l] (solid line). The dashed line represents our 
3 Mj (10 Myr old age) observational constraint and assumes that 
the planet's semimajor axis = its projected separation at the epoch 
of our observations. The dash-dot line is the same, except showing 
our 30 Mj (1 Gyr old age) observational constraint. If HD 15115 
is old, then the object creating the gap must reside between ~ 34- 
45 AU. If the system is young, the allowed semimajor axis range 
shrinks to 40-45 AU. 

would then be constrained to be between 0.87-1" (40-45 
AU). If the system is old, the object must be less than 
~ 30 Mj and its semimajor axis would be constrained to 
0.75-1" (34-45 AU). 

It is certainly possible that we did not detect the pu- 
tative planet because it is currently in front of or behind 
the star. It is also possible that the gap has not been cre- 
ated by a planet, but instead arises from other dynamical 
interactions. Additional high-contrast, high angular res- 
olution imaging of the debris disk would help constrain 
the existence of both the gap and any possible planets. 

5. SUMMARY 

We have presented several intriguing results on the de- 
bris disk surrounding HD 15115. These results are: 

1. We detect east-west asymmetry in the disk mor- 
phology at Ks band, with the western side of 
the disk being a ~ magnitude/arcsecond 2 brighter 
than the eastern side at the same stellocentric dis- 
tances. The asymmetry and brightness differences 
are consistent with results at shorter wavelengths 
and lend additional evidence to the interpretation 
that ISM interactions a re affecting the disk struc- 
ture (jDebeset alj|2009t ). 

2. At L', we detect symmetric and ~ equally bright 
disk structure. 

3. We detect an overall bow- like shape to the disk at 
both Ks band and L' , and the disk appears offset 
from the star to the north by a few AU. We are 
able to explain these observed effects using a model 
disk with a near edge-on inclination, a gap between 
0- l'.'l, large (3 /im) grains, and mostly forward- 
scattering grains (g = 0.5). 

4. The FWHM of the disk at Ks band and at 
V is consistent with the disk FWHM at 1.1 



5. The disk SB profile at Ks band shows evidence for 
a 2er reduction interior to l'.'l (50 AU). Combined 
with SED analysis, this may be a sign that the disk 
has a gap interior to 1" (45 AU). Additional high- 
contrast observations at NIR wavelengths with bet- 
ter inner working angles would help confirm or dis- 
prove the proposition that the disk has a gap. 

6. The Ks - L' disk color is mostly grey for both 
sides of the disk between 1.1-1'.'45 (50-66 AU). 1-10 
fiia grain sizes are the best match to the western 
disk color. 3-10 fiia grain sizes are the best match 
for the eastern side. Given the system's expected 
grain blowout size of 1-3 /um, our dust grain size 
constraint s may support the ISM interaction inter- 
pretation (|Debes et al.l 120091 ). which predicts that 
small grains would be blown to the western side 
of the disk, leaving large, unaffected grains on the 
eastern side. SED analysis, combined with our ob- 
servational evidence for a gap near 1" (45 AU), also 
predicts a dust grain size of ~ 3 /zm. 

7. We do not detect any > 5cr point-sources at V in- 
dicative of planets. We rule out companions more 
massive than ~ 3 Mj beyond 0'.'6 if the star is 
young and more massive than 30 Mj beyond 0'.'7 
if the star is old. Independently we constrain the 
mass of a single, coplanar, low-eccentricity planet 
creating the gap in the disk to be ~ a few M j if it is 
close to the disk edge, and to be in the brown dwarf 
regime if its orbit is closer to the star. Assuming 
the object's semimajor axis = its projected sepa- 
ration at the epoch of our observations, the planet 
would be < 3 Mj and orbit between 0.87-1" (40-45 
AU) if the star is young; if the system is old, the 
object would be less massive than ~ 30 Mj and 
orbit between 0.75-1" (34-45 AU). 
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Fig. 11. — Expected and observed artificial disks at Ks band and L' . Model disks are always vertical. Top left: Ks band image of the 
artificial disk, showing what the disk should look like if unaltered by the LOCI algorithm reduction. Top right: Ks band image of the 
recovered model disk. Bottom left: the same as (a) except at U . Bottom right: the same as (b), except at L' . In all cases, the model disks 
are recovered, though some self-subtraction is evident. 



APPENDIX 
ACCOUNTING FOR DISK SELF-SUBTRACTION BY LOCI 

To account for disk self-subtraction by the LOCI algorithm, at each wavelength we inserted an artificial model disk 
into the raw data, at a PA nearly perpendicular to the known disk, ran the data through the pipeline, and recovered 
the model disks. The artificial disks were set to SB levels comparable t o the real disks at each wavelength, and the 
widths of the artificial and real disks were compar able. Fig. I Hal and Fig. Illbl show the expected and recovered model 
disks at Ks band, respectively; Fig. Illcl and Fig. Illdl show the same at V . 

After recovering the model disks, we compared the observed PA, FWHM, and SB values as a function of distance 
from the star with the expected PA, FWHM, and SB values. The calculations were performed with identical methods 
to the real disk data analysis. At both wavelengths, the expected and observed PA values are consistent with each 
other, therefore no correction was needed (see Fig. I12al and Fig. Il2bj) . At Ks band, an addition to the observed 
FWHM of ~ 0'.'03 was needed to correct the a pparent constant offset (see Fig. I12cp . At L' , a correction was also 
needed, with a value of ~ O'.'ll (see Fig. I12dl) . Both FWHM correction offsets have been included in the FWHM 
analysis of the real disk images. 

To correct for the reduction in disk SB at both wavelengths relative to the expected values, we mapped out the 
self-subtraction as a function of distance and multiplied this into each recove red d isk. The corrected SB values, along 
with the expected and observed, are shown in Fig. I12el (Ks band) and Fig. Il2fl (L 1 ). The SB corrections have been 
included in the SB analysis of the real disk images. 
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Fig. 12. — Top left: Expected (blue) and observed (red) model disk PA vs. distance from the star at Ks band. No correction is needed. 
Top right: the same, except at L . No correction is needed. Middel left: Expected (blue), observed (red), and corrected (green) model disk 
FWHM vs. distance from the star at Ks band. An offset correction of ~ 0'.'03 is needed. Middle right: the same, except at V . An offset 
correction of ~ O'.'ll is needed. Bottom left: Expected (blue), observed (red), and corrected (green) model disk SB vs. distance from the 
star. Bottom right: the same, except at L'. 
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